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ABSTRACT 


Watkins, William Allen. Ph.D., Purdue University, 
Peeenver  lo/Za Awotudy of Lateral Diffusion in Thin Film 
Bouples of Gold and Silver. Major Professor: 
samuel J. Hruska. 

An experimental deposition technique was developed to 
Teodclead thin film couple after the standard bulk diffusion 
Sandwich couple. The initial region of mixing between the 
two sides of the thin film couple was less than 2 microns. 

An empirical conversion parameter, following the bulk 
procedures of Ziebold and Ogilvie, was determined for thin 
films of gold and silver between 500 and 2,000 A. aay 
intensity data from the electron microprobe were converted, 
through use of the parameter, to units of mass concentration. 
Pure and alloy films of gold and silver exhibited linearity 
between the microprobe x-ray intensity and film thickness 
up, bo. 2,000 A, The maximum error in alloy film composition 
using the conversion parameter was shown to be 10 atomic % 
dtr o00 A. 

Thin film couples of pure silver and gold annealed above 
500°C, and at 416°C for times greater than 72 hours, exhibited 
dewetting on the pure gold side. The dewetting of the gold 


film was shown to be in accord with Mullin's theory of grain 


boundary grooving. For Ag-Au couples annealed below 500°C 








xi 
OEZoNewolspOroOsity Occurred, increasing with time and tem- 
pPomerurcoesotednnedal, leading to film scparation near the pure 
Pivemeside OL the Original region of mixing. The porosity 
formation was similar to that reported by investigators of 
bulk sandwich couples of pure Au and pure Ag. 

Incremental thin film couples with initial concentration 
differences of 20 atom% did not exhibit the zone of porosity 
which again was in agreement with investigators using bulk 
couples with small concentration differences. The diffused 
concentration profiles were compared with numerical solutions 
of the initial profile. The numerical technique was based 
On an approximate solution to Fick's second law for one 
dimensional diffusion. 

A series of apparent diffusion coefficients was deter- 
mined for temperatures between 200 and 416°C. A plot of 
ineeversus 1/T was linear with a coefficient of correlation 


Peeve 72.) ihe data fit the following relationship 
Nel 9exet0 9 “exp(-10,900/RT) cm*sec™? 


An apparent diffusion coefficient was also determined 
for a series of pure Ag and Au couples annealed at 416°C and 
compared favorably with the incremental couples' diffusion 
Geetficient at the same temperature. 

Two models were considered in order to attribute the 
eeercmteaittusivity to a single transport mechanism. The 
surface diffusion model was refuted because of the low 


reported value of activation energy (10.9 kcal/mole), in 





X11i1 


comparison with what was theoretically expected and experi- 
mentally reported by other investigators. <A model attribut- 
ing the apparent diffusivity to transport along the film's 
grain boundaries was considered to be likely, if the grain 


boundaries could be assumed in a high state of disorder. 





INTRODUCTION 


Titers tonein thage(< 5,000 A) vacuum deposited metallic 
films is a phenomenon which has received significant atten- 
mEoneonty aurine the last ten years. This is largely due 
Peele wiNGreasine use Of thin films in electronics either 
aoeecomorete thin film Circuits, or more commonly as inter- 
connections and contacts in integrated circuits. 

Ui@MaeenewlasSt &WO years, the Majoritye of thin film 
diffusion studies had been done with multilayer films. 
Multilayer films are produced by the successive deposition 
Of material stacked layer by layer, one on top of the other, 
as shown in Figure 1. In a layered film, deposited in high 
vacuum (< 10 °torr), each interface region will be nearly 
perfect. For example, there will be atomic contact between 
any two superimposed deposits. Diffusion occurs in one 
direction only, and the limiting boundaries of the system 
Memo nvetaretlaneconthne direction of the highest concentra- 
tion gradient. Several investigators have devised experi- 
mental methods to detect the subsequent diffusion over the 
extremely short distances relative to that of bulk diffusion. 
Two comprehensive and fully referenced review articles have 
been written by C. Weaver [1, 2] and bring the reader up 


through 1970 in all work related to diffusion in thin films. 





METAL — METAL 
INTERFACES 


DIFFUSION DIRECTIONS 


EEESSSEL. 








Jee ure al 
Multilayer Film Geometry Using Alternating 


Films of Metal A and Metal B 





Do 





Fiee~eGomlections and Contacts in inteprated circuitry 
made of dissimilar materials have a possible mode of atomic 
transport which cannot be modeled using the multilayer film 
technique. <A thin film (< 2,000 A) couple made of two films 
Side by side on a suitable substrate has an increase in the 
surface to volume ratio of approximately 10° relative to a 
Peancdnd Size bulk diffustom couple. Atomic transport 
feos cie t1imestiniace-Vapor interface could significantly 
Mientpute coOethe rate of homogenization, especially at 
temperatures below 3 the melting point of the pure components 
Teen SeOuUDTC ii substrate interactions and surface 
morphology changes for various annealing times and tempera- 
ures Can also be studied. 

Patera adaitision. arrtusion along the length of a 
binary thin film couple parallel to the substrate, has not 
had the same degree of attention as that of diffusion normal 
to the substrate (multilayer films). Only four investigators 
[3, 4, 5 and 6] have developed experimental techniques to 
investigate lateral diffusion. Three authors [3, 4 and 5] 
used the overlapping film technique. This method consists 
of alternating film deposition onto a suitable substrate 
intervened by a mask shift parallel to the substrate. The 
resultant film-substrate geometry used by all four investi- 
gators is shown in Figure 2. 

Metermitime deposition, each investigator gave the 


binary thin film couples an anneal in vacuo. The overlap 


a ~~ ee 
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REGION OF 
MIXING 





Film Geometry for Previous Lateral Diffusion Studies; 


A. Monch [3] B. Pines et al. [4] and Johnson and Brown [5] 


C. Smallen [6] 








Pete weOtpy siege b. In Figure 2, homogenizes rapidly duc to 
Siem sior teri 1usion distance perpendicular to the substrate. 
imG@morethce four investigators [3, 4] overlapped films of 
eevee he wandsoc, respectively. Another author [5] studied 
diffusion in both sets of binary systems. Both Ag-Te and 
Ag-Se have an intermetallic AgoX3 reported [7] in studies 
of their equilibrium phase system. All three authors [3, 4, 
Si) reported rapid formation of the intermetallic compound in 
the overlap region. Further, each of these investigators 
studied the subsequent lateral growth of the intermetallic 
Pacer wnonMoeenizdation Of the overlap region. High rates of 
diffusion with low activation energies, relative to bulk 
Meiustonevalues, were reported. [3, 4] attributed these 
hesultsS to transport along grain boundaries, and [5] to 
(imicmonealone the film sSurface-vapor interface. In all 
these three independent studics, the experimental results 
were strongly dependent on the subscquent film geometry. 
ferwatdenssim either of both of the film's thicknesses 
Gauscd variations of two orders of magnitude in the reported 
fet risieneratesss In (5), the Ag had to be deposited first 
Snethe substrate in order for lateral diffusion to occur at 
all. 

Smallen's work [6] in the gold-copper system provided 
the major impetus for the present study. Using the method 
of simultaneous deposition of the gold and copper onto a 


Suitable substrate, he reported rates of diffusion both 





parallel and perpendicular to the substrate. Smallen 
feportcamrares Ol Gdritusion in the lateral direction almost 
seven orders of maynitude faster than in the perpendicular 
direction. An associated activation energy was reported and 
was almost 3} that of the previously reported bulk value. 

A study of the Bi-Te system, which also has a reported 
[8] intermetallic compound in its equilibrium phase systen, 
Was undertaken by this author. In order to more closely 
approximate the sandwich couple, so as to compare more 
favorably with the initial conditions of the bulk, it was 
decided to try to improve the experimental methods of [3, 4, 
Sand 6 | . 

An experimental method was devised and reported [9] in 
order to reduce the overlap region of [3, 4, 5 and 6] from 
the order of 10°? meters, to 10°° meters, and minimize the 
Mibomlcss Variation in the lateral direction. The new 
fnetiyetOL eas tnin tilm diffusion couple used in the 
Prono sewohk 1s shown in Figure 3. Note that the overlap 
region of Figure 2 is now called a region of mixing, which 
more closely represents the physics of the new experimental 
deposition process [9]. 

The formation of an intermetallic compound in the 
region of mixing was observed in couples made of pure 
bismuth and tellurium. However, after formation of the 
intermetallic compound, no lateral diffusion was observed. 


A summary of the work accomplished on this system, with 
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DIFFUSION 
DIRECTIONS 
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Geometry for Thin Film Couples Used in Present Study 





eee wemimenidl results, and a discussion of the same, is 
fOouna in Appendix A. 

A second study was undertaken in the Ag-Au system, 
Mitel seemtirely miscible throughout the entire composition 
range. In this study, the region of mixing was further 
reduced to the order of 10°° meters as the experimental 
procedure was further improved. (See Experimental Procedure 
pection. } 

The initial phase of the second study involved couples 
of pure Au and Ag. During the anneals surface morphology 
changes occurred in the couple. The changes were similar 
bo those previously reported in studies of bulk Ag-Au 
Pouples, Nowever they occurred to a greater extent and 
interrupted the continuity of the couple. Because of the 
film separation, quantitative rates of diffusion could not 
be determined as anticipated. However, similarities between 
the concentration profiles for various annealing times led 
to an alternate technique; this allowed values of diffusivity 
Pompe determined. 

Johnson [38], using bulk Ag-Au couples with small 
initial concentration differences, did not observe these 
morphology changes. Accordingly, an experimental procedure 
was devised to make thin couples with an initial concentra- 


0 


4. As with Johnson's results, 


nm 


BIOneamtrTerence of 20 atomic 


no morphology changes were observed in the present study. 








A numerical technique was devised to treat experi- 
(i Wewevecetermined initial concentration profiles. The 
numerical solutions were compared with the experimental 
concentration profiles after various times and temperatures 
of anneal. Values of chemical diffusivity were determined 
as a function of annealing temperature. From the comparison 
Seecne numerical solutions with the diffused concentration 
mrotiles, an e€quation for the apparent diffusivity for the 
couple was determined. Possible atomic transport mechanisms 
were investigated over the working range of the present 
study. 

Lastly, a calibration study, involving the fabrication 
and use of Ag-Au thin film standards was undertaken. 

kinpinredal working curves over the entire composition 
range were determined so that electron micro-analysis could 
be directly applied to deduce mass concentration variations 


iimene dittlision studies. 
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BACKGROUND AND LITERATURE SURVEY 


Two separate investigations into the literature had 
to be undertaken in order to carry out the experimental 
work of this author. 

The first investigation focused on a method to analyze 
the annealed thin film couple and from this, deduce mass 
concentration variations as a function of lateral distance, 
time, and temperature of anneal. 

The second involved an investigation into the previous 
diffusion studies in the Ag-Au system in order to compare 


Poetieetne C€xperimental results of the present thin film study. 


Empirical Method for Thin Film Analysis 

The experimental method of [3, 4, 5] involved the 
optical observation of the growth of the diffusion zone as 
a function of time and temperature. In all three cases, the 
diffusion zone was an intermetallic compound with a signifi- 
cant reflectivity Weererence relative to the pure components. 
Smallen [6], however, worked in the gold-copper system which 
is completely miscible throughout the composition range. 
In order to determine atomic concentration profiles in his 
couples he used the electron microprobe. By scanning his 


samples, under the electron beam, he generated x-ray 
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Pitemstityvedata sds a funetion of lateral distance. In order 
MemcOlvenrtUstiic X-ray Spectrometric duta to atomic concen- 
Moron, ne macae the wmather crude assumption that the 
miative xX -1ay Intensmuey,, corrected for background, was 
proportional to atomic concentration. 

Due to the similar work in miscible systems by Smallen 
[6], and to the present work of this author, the electron 
microprobe was used to determine mass concentration changes 
as a function of lateral distance, time, and temperature of 
Mmuled! —wiOWeVer, an InveStigation was necessary to deter- 
mine the relationships between the x-ray spectrometric data 
generated from a thin film sample, under the electron bean, 
and the specimen's mass concentration. 

The important variables relating the emitting x-ray 
imicensity tO the chemical composition of the beam target, 
were first reported by Castaing [17]. The full correction 


for the measured relative intensity in the bulk is 


where {f£(x)/f°(x)] is the absorption correction factor 


gine 9K MomrEne COrnecti1on for secondary 


¢] 
fluorescence 


ES Womthem poe kscatter-1Onization Correction 


is the weight fraction A in A-B alloy 





: 


az 


Nikos is the measured x-ray intensity of the 
Characteristic line for element A 
(corrected for background) of the A-B 
BEL Ory. 

ee a is the measured x-ray intensity of the 


characteristic line for element A 
(corrected for background) of a sample 
of pure A. 
Ziebold and Ogilvie [18] made a series of experimental 
observations between Ky and Cy mob wdifierene binary 
Systems, including the Ag-Au system. They found that their 
experimental data in all 17 systems displayed linearity 
qa» Versus Cy. The analyti- 


cal straight line function is written 


when plotted in the form of C,/K 


Cy/Ky = (1 - Aanp)Cy- (2) 


The single parameter needed to describe the empirical work- 
ing curves of [18] is the limiting value of Cy/Ky, when Cy 
Peesmtoezeno. the subscript, AB, indicates that the con- 
version parameter applies to the measurement of characteris- 
tic x-rays of element A in binary A-B alloys. 

Rearrangement of equation (2) gives a more symmetrical 


Owen 2 


ee 2 ae (3) 


i? 
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Figures 4 and 5 show the experimentally determined bulk 
relationships of [18] between K and C in the silver-gold 
system. A hand calculation of the values of the three 
correction parameters of equation (1), as a function of 
composition, was made using semi-theoretical formulae, 
based on the three correction factors of equation (1). 
Pipoperating voltage of 25 kv of the electron microprobe, 
with a 52.5° spectrometer takeoff angle, (parameters used 
in the present experimental work) was assumed in the calcu- 
Pattomes ime rigure 6, the results of the hand calculation 
are plotted in terms of the variables of equation (3). 

In the case of thin film analysis, if the linear 
analytical function of equation (2) is assumed to hold 
POretnin films inthe Ag-Au system, then the empirical con- 
version parameters (an sau? a nuAg? must be determined 
separately from the previously reported bulk values. 

Hutchins [19] showed that the thin film x-ray intensity, 
normalized with respect to a bulk sample of the film material, 
Sameber expressed aS three factors: 


I 


FS 1 
le) Soe Cl nT 62) (4) 
bulk aes Learn exc. . 


where Inc = the measured intensity for the film on a 


particular substrate 


Lake = the measured intensity for an infinitely thick 
film (with respect to the electron beam) 
Se-esupstrate correction 
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INTENSITY 





MEASURED RELATIVE 


O 0.5 1.0 
WEIGHT FRACTION GOLD 


Figure 4 
Experimental Microprobe Calibration Data for the Bulk 
Silver-Gold System at 30 kv and 15.5° Spectrometer 
Takeoff Angle. 


Zeibold and Ogilvie [18] 








Ca/Ka 





WEIGHT —_—"" GOLD 
Ca 


Figure 5 
Functions Cy/Ky for Bulk Silver-CGold System. 
Micromropc at 50 ky and 15.5° Spectrometer 
Takeoff Angle. 


Zeibold and Ogilvie [18] 


15 
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ee, 
Cr 


Fe abe d 5e(S: 
Results of Hand Calculation Using Semi-Theoretical Formulae, 
Based on Equation (1) Plotted in Form of Equation (3). 


Microprobe at 25 kv with 52.5° Spectrometer Takeoff Angle. 





Wye 


I, = the measured intensity for the film without a 


substrate 


eee Te Aes miemexecired Intensity for an infinitely 
eee Kea lm 
f(x) = the x-ray absorption correction for film 


thickness, t. 
For a particular set of operating conditions of the 
electron probe and film thickness, t, the measured relative 


intensity can be written as 


I 
a (Ss 
pure A(FS) 


where I alloy (FS) the measured intensity for an A-B alloy 
film on a particular substrate at a film 
thickness t 
H beae ACES) = the measured intensity for a film of 
pUyeeowONea particular substrate at a 
Pein threkness  t. 
Solving equation (4) for Inc, and substituting into the 
numerator and denominator of equation (5) with appropriate 
subscripts 


I 
K, (film thickness t) = EL) 


Ture A(F) 


Thure A bulk exc. 1%)p ot tate (6) 


I 
I alloy bulk exc. £(X) y bulk pure A 
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Ziebold and Ogilvie [18] were unable to justify their 
experimental observations leading to equation (2) on the 
mieonetlecalegrounds of the bulk relationship of equation (1). 
This author will likewise make no claim to the prediction of 
Pelinear relationship in the present thin film data plotted 
in the form of equation (2), based on equation (6). A 
Straight line relationship, within the limits of the repro- 
ducibility of the experimental data, between the variables 
of equations (2) or (3) will justify the use of the empirical 
conversion parameter to deduce weight fraction composition 
from the measured relative intensity of the electron micro- 
probe. 

To convert from weight fraction composition (C,) to 
mass concentration involves knowledge of the density varia- 
tion of gold-silver alloys across the entire composition 
range. 

In Figure 7, the Vegard's law deviation is shown for 
the silver-gold system. Using the lattice parameter data of 
Sachs and Werts [20], and the relationship derived in 
Barrett and Massalski [22] for the density in a substitutional 


alloy system, 


o(AgAu) = nA/VN (7) 


where N Avagadro's number 
V = volume of the unit cell 


N = number of atoms/unit cell 


>| 
i 


mean atomic weight, 
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Picmcenocityeat poOld-silvyer alloys as a function of atomic 


concentration was determined and is shown in Figure 8. 


PONE 0S 0 


At the present time there has been no reported study of 
thin film diffusion in the gold-silver system. 

The initial studies of chemical bulk diffusion in the 
silver-gold system was undertaken by [23, 24 and 25]. In 
all three works, the method of chemical etch limit was used 
to determine average diffusion coefficients over the entire 
composition range. Johnson [38], who doubted the validity 
of the chemical etch limit method, studied chemical diffusion 
in the silver-gold system using sandwich couples with initial 
concentration differences of less than 20 atomic %. Concen- 
tration-penetration profiles were determined by the now 


accepted method of sectioning. Using Fick's Second Law, 


Can? dC 


Johnson assumed that, due to the use of couples with small 
initial-concentration gradients, D was not a function of 
composition over the working range of his couples. Therefore 


he wrote 


Q 
- 
@ 
ied 
2. 


(9) 


| 





Q> 
Ge 
Q 
vad 
No 


and solved equation (8) using the standard initial and 


boundary conditions, 
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Variation of Alloy Density for Bulk Alloys in Gold-Silver 


System Using Relationship Derived in Barret and Massalski [22] 
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Johnson compared his experimental profiles with those of 
equation (10) to determine D, and also plotted his results 
on probability paper to justify his assumption of D # £f(C). 
Figure 9 shows Johnson's results for an alloy containing 
49.2 atomic % gold, along with the results of [235, 24 and 
25]. Jost's [25] data, taken at lower temperatures, were 
extrapolated by Johnson [38] up to his experimental tempera- 
ture range. The data of an investigator named Braune are 
reported in Johnson's work and shown in Figure 9, however, 
he does not reference him in his paper. The straight line 
through Johnson's data can be represented by the following 
equation, 


_ -41,700 2 x 
D. = 0.14 exp [tg] cm“sec °. (11) 


Included in his work were determinations of the gold 
and silver self diffusion coefficients using radioactive 
tracers. Over the same temperature range they are repre- 


sented by the following equations, 
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Dag = 0.39 exp| a7 | Gmosee — Ole 
Dae = (UB exp (te?) em secs (ls 


Pevetyecxtensive Study of self diffusion in silver-gold 
solid solutions was done by Mallard et al. [26], also using 
radioactive tracer sectioning techniques. They reported 
that the silver tracer moved between 2 and 3 times faster 
than the gold tracer, across the entire composition range. 
They also report a variation of less than } of an order of 
magnitude in the diffusivities of both the gold and silver 
tracers across the entire composition range. 

Baluffi and Seigle [27] studied diffusion in bi-metal 
warpor-solid couples of zinc-alpha brass, copper-nickel and 
gold-silver. In their investigations of vapor-solid couples 
@fegeld-silver they observed the formation of 13% volume 
porosity, by lineal analysis, in the diffusion zone after 
93 hours at 940°C. In order to elucidate the cause of the 
porosity formation, they made an additional study of pure 
gold-pure silver sandwich couples. The porosity again 
appeared in the sandwich couple and its formation was 
evident at times as early as 1 hour, increasing with time 
up to 50 volume %. Similar observations on gold-silver 
sandwich couples were made by Seitz and Kotterman [28]. 
Both investigators [27, 28] observed the porosity formation 
on the pure silver side of the couple. Baluffi and Seigle 


[27] made corrections to account for the observed volume 
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porosity to determine diffusivities by the Matano method. 
Seitz and Kotterman did not make corrections and therefore, 
their data are suspect. Baluffi and Seigle took Laue back 
reflection x-ray photograms of the diffusion zone. They 
Susetved multiple spots and attributed this to polygoniza- 
tion of the diffusion zone. 

In his earlier work, Baluffi [29] determined that 
Silver atoms diffuse into gold more rapidly than gold 
diffuses out, causing a volume increase, with expansion in 
all three directions. Using the theoretical calculations 
of Huntington and Seitz [31], he speculated that diffusion 
by vacancy mechanism coupled with the unequal diffusion 
rates, could account for the observed volume increase. 

Barnes [30] used multiple sandwich couples of thin 
interleaved sheets of copper and nickel to accentuate any 
observed volume increase due to diffusion. He reported that 
there was an increase of volume in the diffusion zone as a 
result of atomic diffusion in the copper-nickel and copper- 
zinc systems that he studied, each system having unequal 
rates of the two diffusing species. He also showed that 
this volume increase is associated with the formation in 
the diffusion zone of voids with faceting. These voids 
appear near the original interface and on that side from 
which there is a net loss of atoms (net gain of vacancies). 
He also conducted an experiment which showed the presence of 


strain in the diffusion zone and another experiment to 


26 


SeierpWcemtnis Strain tO ditfusion. He conducted an x-ray 
investigation like that of Baluffi and Seigle [27], and 
again concluded the presence of polygonization in the diffu- 
Sion zone. Barnes [30] speculated that the difference in 
atomic fluxes of the two diffusing species was compensated 
by a vacancy flux in the direction of the smaller flux. 

The generation and annihilation of vacancies on opposite 
sides Of the couple cause strain in the couple. He sug- 
gested that the voids forming as a result of diffusion, 
imply that some of the large number of vacancies involved 
in the diffusion zone condense out of the lattice. He 
postulated generation and annihilation sites such as grain 
boundaries, polygon boundaries and on the edge component of 
aire locations. 

He, however, did not discuss the magnitude of the devia- 
tion of the equilibrium number of vacancies present in the 
Merrusion zone during the anneal, He explained that the 
growth of the void could be explained by diffusion and 
absorption of some of the excess numbers of vacancies on 
void surfaces. 

Seitz [34] postulated that the process of porosity 
formation is precipitation of supersaturated vacancies 
caused by the unequal diffusion current and occurs by 
heterogeneous nucleation. He calculated a value of NOs "ton 
the lifetime number of vacancy jumps in the diffusion zone 


by using the porosity observations of Barnes and an arbitrary 
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Ile wOLwetOnPethe Critical "relative excess vacancy concen- 
tration to form holes. 

Baluffi [32] supported Seitz's postulate of heterogene- 
ous nucleation to form voids from a supersaturated solution. 
Baluffi, however, pointed out that the high value of super- 
Satubation expected by Seitz could not be compatible with 
the experimental results of parabolic Kirkendall marker 
shifts with time. He considered the problem with an 
improved mathematical treatment, and his results indicated 
that the relative excess vacancy concentration for void 
formation to be < 0.01 (1%). 

Barnes and Mazey [33] using sandwich couples of Cu and 
Ni under a hydrostatic pressure reaffirmed Baluffi's [32] 
conclusion of 1% relative excess vacancy concentration 
necessary to form holes, and also showed experimental evi- 
dence that void formation was by heterogeneous nucleation, 
probably foreign particles or inclusions at the weld inter- 
face. A more recent and complete mathematical treatment of 
the vacancy flux contribution during diffusion and Kirkendall 
shift in binary alloys, is given by Manning [37]. 

During an experimental investigation of the Kirkendall 
shift in the gold-silver system, Ruth [35] observed the 
formation of a zone of small protrusions of the pure gold 
side of the couple. He reported that the growth law of 
the width of this zone is a secondary effect of diffusion. 


The protrusions were observed to form in the vicinity of 
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Mibehesectlons With Crystal Surfaces of <112> lines in {110} 
planes. He explained the observations by proposing that a 


: Pero iimeacendistocation, which has a {112] direction, 


2 
acting as a diffusion pipe. The excess silver atoms, 
caused by the unequal diffusion rates, should migrate along 
meee dislocation line to the surface. 

theme are mo reported studies of grain boundary diffu- 
sion in bulk silver-gold alloys. Hoffman and Turnbull [36], 
however, studied lattice and grain boundary self diffusion 
in silver. Over the temperature range of their study (250- 


500°C) they reported the grain boundary diffusion coeffici- 


me times the grain boundary thickness f1t the equation, 


Dep eee = 125 x 107° exp| 20,700) cm?sec™? (14) 


where y = 5 x107°%cm 


hratnesizes for their specimens varied from 21-60 microns. 
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DESIGN AND CONSTRUCTION OF EQUIPMENT 


The successful deposition of thin films in vacuum 
requires a working understanding of vacuum technology. The 
experimental apparatus must be maintained in an environment 
Suitable for the selective transport of mass over a distance, 
uSing materials which maintain this environment during all 
the operating conditions of the experiment. In order to 
give the reader insight into the many criteria which had to 
Bomsatistied in order to successfully create the three thin 
film geometries, and carry out the required experimental 
conditions for the study of atomic diffusion, detailed 
explanation is given in the following sections. 

The discussion of this author's experimental equipment 
is broken into three sections: the vacuum system, the 
experimental apparatus in vacuo, and the associated external 


equipment. 


Vacuum System 


The environment for the generation of both binary thin 
film couples and thin film standards was a vacuum chamber 
capable of maintaining residual gas pressures of 10° ’torr. 
Primary pumping was obtained using a 4 inch oil diffusion 


pump in conjunction with a mechanical roughing pump. 
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Auxiliary pumping, located in a stainless steel base 
separated from the diffusion pump by a high vacuum valve, 
included a Varian tri-filament titanium sublimation pump 
and a liquid nitrogen cooled cyrogenic baffle. Residual 
gas pressure was measured using a Granville-Phillips dual 
tungsten filament ionization gauge tube with an associated 
Pemres 260 controller. 

A 12" Varian, all stainless steel, spool piece was 
placed on top of the base, using viton gaskets to maintain 
vacuum integrity. Integral to the spool piece are ten 22" 
Conflat ® ports for insertion of necessary feedthroughs to 
maintain electrical and mechanical connectivity between the 
evacuated apparatus and the external test and monitoring 
component, Am 18" high pyrex glass bell jar, also sealed 
with a viton gasket, was placed on top of the spool piece 
completing the vacuum system described above. A sketch of 
the vacuum system with its important components labeled is 
Shown in Figure 10. Photographs of the exterior and 
interior of the vacuum system are displayed in Figures ll 


and 12 respectively. 


Experimental Apparatus in vacuo 


The substrate holder, two sources, andassociated monitor- 
ing and shielding mechanisms were constructed entirely within 
the 12'' diameter spool piece. OFHC copper, type 303 and 304 
stainless steel, aluminum, and 96% polycrystalline dense 


alumina, were used exclusively in the construction of the 
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Schematic of the Vacuum System Used in the Present Study 
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Figure ll 


Photograph of the Vacuum System and Associated 


External Equipment Used in the Present Study 
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Femiee 2 


Photograph of the Experimental Apparatus in vacuo 
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experimental apparatus. These materials have reported [10] 
low (< 10°%torr) vapor pressures at the operating tempera- 
[ices Of the vacuum system (-192°C to 600°C). In addition, 
wherever possible, leak paths were provided for all fastened 
assemblies to minimize the outgassing of virtual leaks. 

The substrate holder, a sketch of which is shown in 
Figure 13, was made from a 3" x 2" x $" block of OFHC copper. 
Both large faces of the copper block were milled to within 
SOM ee lwo similarly milled copper restraining plates were 
used to hold the substrate uniformly against the copper 
block. All materials, except the heater and insulating tubes, 
which will be described later, were constructed out of OFHC 
@epoer. tlhe use of copper exclusively on the substrate 
holder eliminated differential stress buildup during heating 
of the block, maintaining the critical alignment between the 
various support mechanisms both on and off the substrate 
holder. 

A copper quench tube, making two passes through the 
substrate holder, was connected to two stainless steel Swage- 
lock fittings. The fittings completed the quench line, 
through attachment to two stainless steel quench gas delivery 
tubes connected externally to the gas source by a stainless 
steel dual tube feedthrough. In addition, the fittings 
allowed removal of the substrate holder, by sealing together 
through the use of replaceable copper gaskets. The razor 


blade support mechanism was made of two sets of two small 
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Schematic of Bottom View of the Substrate Holder 


with Substrate and Razor Blade Mounted in Place 
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copper blocks, squared true to the substrate holder's milled 
face. In each set where the faces held the razor blade, a 
.002'"' reset was milled on each face to accommodate the razor 
blade thickness. Copper screws held together the two sets 
of blocks in each set, and as a unit to the substrate holder. 

Chromel-alumel thermocouples used to record the sub- 
Strate temperature were press fit by screw heads onto the 
surface of the substrate away from the line of sight of the 
two sources. 

Four (4) tubes, symmetrically placed above and through 
Mitcesuibstrate and copper block respectively, were filled 
with sleeves of 99.7 wt. % alumina of 96% polycrystalline 
density. The sleeves provideddelectrical insulation between 
micwcoppen DIOCKk and the .010'’ coiled tungsten wire, which 
was interwoven between the four sleeves, providing a heat 
source to the substrate. Alumina of lesser 
density was found to absorb water vapor readily at atmos- 
pheric pressure and subsequently outgas heavily during 
heatup in vacuum. 

The deposition of gold and silver was carried out in 
two independently controlled resistance heated boats 
purchased from the R. D. Mathis Company. A .010"" molybdenum 
boat with a .010" coating of alumina was used for the gold 
deposition, and a .010'' tantalum boat was used for the silver. 
The tungsten boat more commonly used for gold deposition in 


flash evaporation systems was found to be partially wet by 
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Mico givingmerratic rates of deposition, and possible 
Pericinrcatton of the deposited film. 

Figure 14 is a sketch of the approximate orientation 
relationships between the sources, substrate holder, shutter 
mechanism, and crystal mass deposition monitors. 

The shutter mechanism moves, according to the plane of 
the sketch in and out of the paper, on two aluminum tracks 
fastened to the spool piece. The shutter was designed with 
the winged contour to provide, simultaneously, controlled 
interruption of the gold and silver vapor to the substrate, 
continuous monitoring of each vapor flux by the respective 
crystal monitors, and vibrationless transit with respect to 
the substrate-razor blade-flux barrier alignment. 

As indicated in Figure 14, the flux barrier between the 
twoesources 1s completed by a stainless steel sheet. In 
order to compensate for expansion of the substrate holder 
during heating, and adjustment of the edge of the sheet true 
to the edge of the razor blade, the flux barrier is adjust- 
lebeavetmmcaliyv and horizontally, respectively. 

The radiant transfer of heat to the vacuum chamber 
Wels sendesurrollding apparatus during heating of the sub- 
strate holder assembly was minimized by the construction 
of a two piece (for disassembly) copper box surrounding the 
substrate holder on three sides (open side facing sources). 
The addition of the copper shroud served a dual purpose. 


It maintained the internal surfaces of the vacuum surfaces 
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Figure 14 
Sketch of the Approximate Orientation Relationships 
between the Sources, Substrate Holder, Shutter Mechanism 


and Crystal Mass Deposition Monitors 
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below 100°C, preventing outgassing of the same reducing the 
low pressure capability of the system. Secondly, it assisted 
in maintaining an isothermal condition for the substrate. 
Two Sloan water cooled ultra-high vacuum remote crystal 
holders were rigidly mounted in line of sight, one to each 
source. Five (5) MC gold and silver coated quartz, AT cut 
round crystals were used in the holder for monitoring the 
gold and silver source respectively. Each crystal gives 
thickness monitoring up to 60,000 A of aluminum or mass 
equivalent. The round crystals used were found to be 
Superior to the square crystals used in early experimenta- 
tion. The uniform stress pattern generated in the round 
crystal during mass deposition eliminated false (unrelated 
to mass deposited) frequency shifts common to the square 
type. A radiant heat shield placed around the crystal 
minimized false frequency shifts of the crystals due to 


induced thermal stress. 


External Equipment 


Automatic deposition control was achieved using two 
Sloan Omni II Deposition Control Monitors. While each 
source was controlled to give the preset deposition rate, 
the rate and thickness output was continuously monitored 
by 4 Rustrak chart recorders. 

Substrate temperature was monitored by a Leeds and 


Northrup millivolt potentiometer, Sargent Model SRG chart 
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recorder and a Leeds and Northrup M Line Model C-l, C.A.T. 
controller. 

[meconnuncuron wien the C.A.1. controller, a Leeds and 
Northrup M Line Model C-2 low level set point unit, silicon 
controlled rectifier, and tungsten heating element in the 
substrate holder assembly was used to control the substrate 
temperature during the anneals to within + 1°C. 

A variac auto transformer fed by 115 VAC from a con- 
stant voltage transformer was used for initial heatup and 
outgassing of the substrate and was switched in or out of 
Miemieating Circuié with a DPDT high current switch. 

The bottled helium quench gas was of commercial purity, 
regulated, and provided to the appropriate feedthrough, 


maintaining the substrate holder assembly, by tygon tubing. 
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SUBSTRATE 


When studying the physical phenomena in thin films, 
the substrate or restraining medium on which the films are 
created is of signal importance. The substrate in many 
Cases may determine whether or not a particular physical 
phenomena in the film or films occurs. It was with this 
in mind that this author has written this section, and 
mvaced it into two parts, the selection of the substrate, 
and the subsequent substrate preparation prior to actual 


film deposition. 


Substrate Selection 

Single crystal silicon wafers grown in the [111] from 
the melt were purchased from the General Diode Corporation. 
The wafers were 1.5" in diameter and 10-15 mills in thick- 
ness. As received, one side was mechanically polished, the 
Cemer etched. 

The selection of the silicon wafer as a substrate for 
all experimental work in the silver-gold system was based 
on the following criteria: Firstly, the substrate must 
have a surface flatness and roughness so to consistently 
reproduce the desired film geometry with respect to the 


region of mixing and uniform thickness. Deviations from 
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Efiatness in the substrate leads to variations in the distance 
between the razor blade edge and substrate. This coupled 
with the intrknsic variations in the razor blade edge height 
eameltead £6 umbra effects resulting in formation of a gap 

in lieu of a region of mixing between the two films. The 
surface roughness of a substrate leads to surface protrusions 
of high frequency. The subsequent decay by surface self 
diffusion is considered an undesirable variable for the 
present experimental method. Doctors Anderson and Neudeck 
{11] compared nine common substrate materials and found the 
mechanically polished silicon wafer to exhibit a flatness 

on the order of 0.1 part/thousand, and a surface roughness 
Sumtne Order sor, or below, the resolution limit (= 50 K) Ok 
the measuring instrument [12]. 

Secondly, the thickness of the substrate was to be 
minimized as much as possible to reduce the temperature 
gradient maintained between the substrate holder where the 
temperature was measured and the surface of the substrate 
Supporting the films. A 1.6 mm thick soft glass substrate 
clamped to the smooth face of a copper block held at 90°C, 
exhibited a 20°C temperature difference across the substrate 
thickness [13]. The silicon wafer, with 1,500 A of Si0, 
grown on the surface, is almost an order of magnitude 
thinner than soft microscope glass. A temperature difference 
of 1°C, or less, across the wafer's thickness was experi- 


mentally observed by this author when the substrate holder 
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was held at 416°C, and uniform clamping pressure by the Cu 
restraining plates was applied. 

Thirdly, a substrate with a low average atomic number 
was desirous for quantitative thin film measurement by 
electron probe microanalysis. High energy electrons back- 
scattered from the substrate and the continuous spectrum of 
the substrate, will enhance the characteristic (Au and Ag), 
X-ray intensity produced by the primary electron beam. Both 
the backscattered electron fraction and continuous spectrum 
intensity increase with atomic number [15]. The enhancement 
of the background intensity limits the theoretical limit of 
detection of trace elements. Campbell [14], using fluores- 
Some x-ray Spectrographic analysis defined the minimum 
detectable limit as "that concentration or amount required 
to give a line intensity above background for counting times 
of 10 minutes." For smaller counting times, since the 
standard deviation is a square root function, the loss in 
sensitivity goes as a square of the factor decrease in count- 
ing time. Birks [16], comparing the methods of Campbell 
with the electron probe, showed that due to the smaller 
counting rates of the probe the limits of detectability are 
2 orders of magnitude higher (1 ppm vs. 100 ppm). 

An additional consideration for electron microanalysis 
of thin films is the presence of materials in the substrate 
which cause fluorescence excitation by line radiation. Tin 


and thorium generate a high fluorescent yield for Ag and Au 
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mespectively > however these materials are not found as major 
constituents in most common substrate materials. 

Fourthly, the linear coefficient of expansion of the 
substrate must be as close as possible to those of Ag and 
Peco eprevent Shivering Or crazing to develop in the filn, 
during the quench to room temperature at the conclusion of 
the anneal. Silicon, while not an ideal match, with a 
mie ar coeitiicient of expansion of 7.3 microinch/inch’C, 
Wemees 1422 and 19.7 for the gold and silver respectively, 
gave no experimental observations of shivering or crazing 
mom quenches up to 500°C. 

The above criteria are all satisfied by the selection 
of a silicon crystal wafer as a substrate. Individual 
Criteria may be met to a greater degree by a different type 
of substrate, however none was found that satisfied all four 
Criteria to the degree as that of the silicon wafer. 

Two criteria were not met by the silicon wafer; however, 
they were satisfied by the growth of a thin film of SiO.” 
by means of oxidation of the wafer. 

The growth of the amorphous layer of Si0, allowed the 
deposition of polycrystalline thin films. Secondly, it acts 
as a diffusion barrier, allowing negligible diffusion perpen- 
dicular to the substrate for the times and temperatures of 
this ‘work. 


* The oxidized layer contains molecules of SiO, Si02, and 
loosely bonded oxygen radicals. 
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Substrate Preparation 


All substrates were cleaned by immersion in hot 
sulfuric and nitric acid for 15 and 10 minutes respectively. 
A rinse in.deionized water was followed by several rinses 
in reagent grade acetone and trichloroethylene. The sub- 
Strates were then placed into a constant temperature zone 
of a tubular furnace fed with dry oxygen at 1,000°C, where 
eo 0.0 A of S10. were grown onto the silicon wafer. From the 
sensitive color change due to selective light absorption of 
enemoxidized layer, its thickness was visually determined to 
be uniform over the wafer to well within 100 A. 

The wafers were then masked and placed into a MRC 
Sputtering chamber, where a 5,000 A thick annulus of tantalum 
Pisiean anner and outer diameter of 1" and 1.5" respectively, 
Meaeesputtered onto the oxidized silicon wafer. The sputter- 


ing parameters were as follows: 


System Base Pressure NO ito rt 

Sputtering Gas Argon 

Sputtering Gas Pressure Sar UN torr 
(pressure controlled) 

Cathode Voltage 4 kv D.C. (regulated) 

Current 90 milliamps 


Deposition Rate 250 A/minute 
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The 5,000 A thick annulus maintained a fixed distance 
between the razor blade edge resting on the tantalum surface 
and the S10, surface. This fixed distance for each experi- 
mental run gave a reproducible region of mixing in the 
resultant film geometry. Figure 15 shows the top and edge 
view of the substrate after the preparation process 


described above. 
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EXPERIMENTAL PROCEDURE 


Three different film deposition methods were used during 
the course of the experimental work of this author. The pro- 
cedures used for the vacuum deposition of thin film standards, 
thin film couples, of pure gold and silver, and an incremental 


couple with an initial concentration difference of 20 atomic 


g 


6, Will be described in three subsections. 


Thin Film Standards 

Oxidized Silicon wafers, with no Ta annulus sputtered 
on the surface, were used as substrates. The 1.5'' wafers 
were cut into small sections using a diamond scribe, and 
mounted onto the substrate holder. The available area of 
deposition was further reduced by masking with a thin .010" 
Ta sheet. This provided a sharp edge for thickness measure- 
ment, and a negligible variation in thickness over the film 
area. The razor blade was not used and the stainless steel 
barrier was lowered away from the substrate holder, its only 
use was to provide a line of sight barrier between the two 
sources to prevent cross contamination. The system was then 
evacuated to 5 x 10°-’torr while the substrate was outgassed 
by heating to 150°C, driving off adsorbed water and other 


volatile impurities on the surface. 
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The sources were heated slowly over a period of 5 
minutes, followed by a 2 minute soak cycle just below 
deposition temperature to bring each source to a uniform 
memperature prior to deposition. 

Deposition was carried out simultaneously, the rates 
(A/sec) of each source adjusted, prior to opening the 
shutter, by the control monitors, so that the time of 
deposition for both sources was equal, and yet still main- 
taining the necessary mass ratio to give the film concentra- 
tion desired. Films of 500, 1,000, 1,500 and 2,000 A in 
thickness, each with compositions of 0, 20, 40, 60 and 100 
atomic 4%, were deposited. A detailed explanation of the 
method of achieving film composition is given in Appendix B. 

Gold and silver of 99.99% purity were used as source 
feeerial £0r deposition of thin film standards. 

Film thickness was measured by a Dektak. Four (4) 
measurements of the edge height were taken on each film, 
with the average thickness reported. 

X-ray emission data were taken with the ARL-AMX Electron 
Microprobe. Intensity measurements were obtained with the 
electron microprobe, operated with an emission current of 
100 pamps, beam current of 0.5 pamps, and an operating 
voltage of 25 kv. The intensities of the gold and silver 
La lines were measured with a LiF and ADP crystal, respec- 
tively. Two (2) measurements for a fixed time count of 10 


seconds were taken at five different locations on each film. 


cd = —- 2s 





apnea ii 





/ 
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Background measurements were taken by moving the spectrome- 


hems Off the peak position. 


Thin-Film Couple (Ag-Au) 

Oxidized Silicon wafers with the Ta annulus sputtered 
onto the surface were used as substrates. The substrates 
were mounted onto the substrate holder and clamped firmly 
metoeplace by the copper restraining plates. A Personna 
74 double edge razor blade, with its commercial stamping 
removed, was cleaned by ultrasonic agitation in reagent 
grade acetone and trichloroethylene. A small strip of 
.009 mm aluminum foil was spot welded onto the razor blade 
to cover the center slot design. The final assembly was 
Cleaned further in solvents and then placed into the razor 
blade support mechanism. The substrate holder was placed 
Meomtie vacuum system. The source barrier was then 
femiis ced norazontally to be true to the razor blade edge, 
and vertically so that as the substrate was heated and 
expansion occurred, the vertical height between the top of 
the barrier and the razor blade edge was as small as 
possible. 

The vacuum system was pumped down to 5 x 10° ’torr while 
the substrate was heated by the variac to the annealing 
temperature. The sources were preheated as described above. 

While the substrate was held above room temperature, 
helium quench gas was slowly passed through the quench tubes. 


The small flow rate was monitored by terminating the gas 
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line in a flask filled with distilled water. This prevented 
oxidation of the inner surfaces of the copper quench tubes, 
which were exposed to 1 atm. pressure. 

When the substrate was at the annealing temperature, 
the Leeds and Northrup proportional counter was switched 
into the heating circuit, controlling the temperature of 
the substrate to within + 1°C. 

Deposition of the gold and silver, while the substrate 
was held at the annealing temperature, was carried out simul- 
taneously. The rates (5 A/sec) of each source were adjusted 
equally by the crystal monitors, prior to opening the 
shutter, so that the time of deposition and film thickness 
of the gold and silver were the same. Silver and gold of 
99.99% and 99.9999% purity, respectively, were used as 
source material for deposition. 

At the conclusion of the film deposition, the anneal 
was begun. The anneal was terminated by removing the flask 
from the terminus of the gas line and then forcing helium 
gas at a pressure of 12 lbs./sq. inch through the quench 
tubes. The cooling curves were exponential in character, 
representing quench rates of approximately 75-100°C per 
minute, down to substrate temperatures of around 60°C. The 
quench rate lowered to approximately 50°C per minute for the 
two low temperature anneals. 

The substrates were then removed from the vacuum system 


EOGeruntner analysis. 
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Film thickness of both the pure gold and pure silver 
portion of the couple and thickness variation across the 
region of mixing were measured by Dektac. 

X-ray emission data were taken with the ARL-AMX 
Electron Microprobe. Intensity measurements were obtained 
with the electron microprobe operating as described above 
in the subsection on thin film standards. Intensity profiles 
as a function of distance were determined by recording the 
characteristic x-ray intensity while the sample was moved 
under the electron beam at a rate of 4 microns/minute. 

Fixed time counts of 10 seconds were taken as the sample 
was scanned. 

Photomicrographs of the film surface were taken on a 
Leitz orthoplan microscope and a JSM-U3 Scanning Electron 
Microscope. The scanning electron microscope was operated 
at 27 kv with a sample current of 1.89 x 10°*° amps. 

Limited x-ray emission data were also taken with the scanning 
electron microscope, taking advantage of the small spot size 
of the electron beam, and will be noted as the data are 
reviewed. 

Grain sizes of the films were determined using the 
transmission mode of a Siemens Elmiskop 101 electron 
microscope operating at 100 kv, with beam currents of 10 


and 60 microamps. 
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Incremental Couple (Au-20% Ag) 

Oxidized Silicon wafers with the Ta annulus sputtered 
onto the surface were again used as substrates. The thin 
film deposition was carried out in two steps. 

Initially the substrate, mounted onto the substrate 
holder, was placed into the vacuum system without the razor 
blade. The substrate was heated to the annealing tempera- 
ture, and then a layer of gold film was deposited across 
the entire substrate surface. After cooling to room tempera- 
ture, the substrate holder was removed from the vacuum sys- 
tem where the razor blade was mounted into the support 
mechanism. After alignment of the razor blade with the 
barrier mechanism, as in the case of the Ag-Au couple, the 
System was again evacuated. The substrate was again heated 
to the annealing temperature, and films of gold and silver 
were deposited over the initial gold underlayer, in the 
manner of the Ag-Au couple described in the previous sub- 
section. 

Film thickness, x-ray emission data, photomicrographs 
and grain size determinations were taken in the same manner 


as described in the previous subsection. 
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RESULTS 


Thin Film Standards 

The determination of the thin film empirical conversion 
parameter between the measured relative intensity of the 
microprobe and thin film composition, involved the fabrica- 
tion of twenty-four (24) different films. Gold-silver films 
of 0, 20, 40, 60, 80, 100 atomic % silver were each made of 
4 different approximate thicknesses, 500, 1,000, 1,500 and 
2,000 A. Each film thickness was measured on the Dektac 
where the accuracy of measurement is + 50 A. Two Plots Or 
measured characteristic x-ray intensity versus measured 
film thickness, were made for each film composition. The 
two plots were of the characteristic x-ray data of the La 
lines taken from the gold and silver spectrometers. A 
linear least squares fit was made on each intensity versus 
tiackness plot. 

From the resultant linear analytical equation for each 
plot, the film x-ray intensity was determined for numerical 
thicknesses of 500, 1,000, 1,500 and 2,000 A. Ky and Ky, 
for each composition and film thickness, in the form of 
equation (1), were determined after correcting both the 


pure and alloy films for their respective backgrounds. 
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The backgrounds of the pure Au and Ag films were deter- 
mined by moving the spectrometers off each side of the 
Characteristic La peak of the Au and Ag line respectively, 
and averaging the two values. The alloy background was 


determined by the following relationship: 


B alloy 7 Boure Au Cau i Boure Ag Cag: be) 
where = aillles = the background of an alloy film of 
thickness, t. 
Bure in the background of a pure A film of 
thickness, t. 
Ge a OMmie traction. 
The values of Rag? Kay and Cag? Cay respectively, were 


plotted in the form of equation (3), where another linear 
least squares fit was made, the slope representing the value 
of the empirical conversion parameter, 4AgAu and anuAg 
Mespectively. the process was repeated for each film 
thickness. In addition, a linear least squares fit was 

Maae stor all the data of each spectrometer to determine 

an average value of the empirical conversion parameter for 
film thicknesses between 500 and 2,000 A. 

The statistical calculations discussed above were 
accomplished by the computer program in Appendix C. Tables 
1 and 2 are the direct computer output printed on the Gould 
Picctrostitte printer, displaying the results of the cali- 
bration data taken off the silver and gold spectrometer 


pespective ly. 
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Incremental Couple (80%Au-Au) 
Table 3 summarizes the experimental runs completed 
during the study of chemical diffusion in the incremental 
couples. The couples consisted of a pure gold film layered 


side by side to an alloy film of equal thickness, and with 


9 


a composition of 80 and 20 atomic % gold and silver, 
mespectively. The film thickness for all the runs was 
approximately 1,500 A; actual measurements made on each 
film using the Dektac are reported in Table 3. 

Five (5) couples were annealed for 24 hours at tempera- 
mures Of 200, 250, 300, 350 and 416°C. Two additional 
couples were annealed at 416°C for 12 and 48 hours. 

One (1) couple (Run 227) was not annealed and represents 
the initial (time = 0) concentration profile. Figure 16 
shows the concentration-distance profile generated from the 
computer plotting program of Appendix D. The computer pro- 
gram reads in the x-ray intensity data from the microprobe 
as a function of lateral distance (Stepsize 0.8 microns), 
calculates the measured relative intensity corrected for 
background, and converts Kas to units of mass concentration 
(grams Ag/cc) through the empirical conversion parameter 
and appropriate value of alloy density. 

Figures 17 through 23 show the concentration-distance 
profiles, generated from the same computer program, of the 
annealed couples. In order to compare the annealed profiles 


with the initial profile, and thus determine values of 
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ible ss. oummary Of Experimental Conditions for Runs 224-232 


Incremental Couples (80Au-Au) 


Fume lémperature Time Film 
No. of Anneal of Anneal Thickness Comments 
aac) * (hrs) (A) 
224 416 24.0 1500+ No interface contact 
225 300 24.0 1480 
270 350 240 S25 
2] 25 0.0 1490 
22% 416 24.0 1480 
229 416 48.0 1440 
230 250 Za, 0 LS ZS 
sl 416 12.0 1530 
LEM 200 245.0 1500 
cae | C 
0 
+ SOA 


f not measured 
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chemical diffusivity as a function of temperature, the 
following analytical procedure was used. 

in the case of the thin film incremental couples, it 
was assumed that Fick's second law for one dimensional 
diffusion is a valid approximation to the measurement of the 
rate of change of concentration with respect to time, and 


1s written; 


ales 
CEG) 
iN 
=, 
Q 
Pa 
Nh 





(1S3 


where D # D(C); D # D(x) 


C = mass concentration 
t = time 
x = distance 


The solution to equation (16), if we have a thin (with 
respect to the lateral direction) strip centered at x = Q, 
With an initial thickness (with respect to the lateral 
merection) 6 and concentration C, is 


mee (0-0 is ne 
C(x,t) =§—:}e exp|- | Cle 
2ViDt 4Dt 


A thin strip centered at x = 6, with an initial thick- 


ness 6 and concentration C2, gives 


Ey CoO 06 _ Cane 18 
Cie + ae exp| oe (<8) 


For an arbitrary collection of thin strips of thickness 6, 


for example, and arbitrary initial profile 
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C(x, t) 7 Cn ; = C_ (nd ,0)6 | (xen )* 
x; = ix, ot). = —_—————- €XD | - =a 19 
n=-o | n=-© 2/nDt : Dt | ca 


The thin film incremental couple was assumed to 
represent a collection of thin strips of lateral thickness 
6, and perpendicular thickness t (measured film thickness). 
Mies ltateral strip thickness 6 was chosen to be 0.8 microns, 
the distance between experimental data points. For the case 
Seeche initial concentration profile of Figure 16, the con- 
centration of each thin strip was determined by the follow- 


ime relationship 
C_ (nd ,0) = Ic, (n8 ,0) + Ca, (041) 6,0)]/2 (20) 


or the average concentration between two successive experi- 
mental concentration data points. 

A computer program was written to numerically solve 
equation (19), for time = 24 hours, and a series of 
expected values of D, using the experimental data points 
Gace initial concentration profile, and is shown in 
Appendix E. 

Viecie Dasis Of the numerical results of the computer 
solution to equation (19) the linear least squares slope of 
tiesmumerically solved concentration profile, between the 
values of 1.56 and 0.56 grams Ag/cc, was chosen as the test 
parameter to compare the experimental profiles of Figures 


ieee Ouplie 2 3. 


_ - (a 
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The linear least squares slope was determined using the 
Seatis tical procedures of Volk {39}, and included a test 
for the standard error of the estimate, sum squares deviation 
ieamy (concentration), coefficient of correction, and variance 
of slope. The statistical procedures are written into both 
the computer plotting program for the experimental profiles, 
(Appendix D), and the computer program which solves equation 
(19) (Appendix E). 

Table 4 summarizes the statistical calculations for the 
annealed profiles and compares them with the numerical solu- 
muons Of Equation (19). 

Frome the comparison of slopes of the annealed profiles 
with the numerical slopes of the solution of equation (19), 
feues of Dl were determined. The concentration profiles 
generated from the numerical solution of Run 227 using equa- 
tion (19) are drawn for the appropriate value of D on each 
esperimenital concentration profile shown in Figures 18-25. 
Figure 24 shows a plot of ln D versus 1/T for the 24 hour 
anneals at 200, 250, 300, 350 and 416°C. The experimental 


data fit the equation 
Doo ee) sills exp|- 10 .f00 cm2sec7! al, 


Simaddition, the values of D, for the 12 and 48 hour 
annealed profiles at 416°C, are reported in Table 4. These 
values were determined in a manner similar to the 24 hour 
anneal, with the slope of their respective concentration 


profiles between 1.56 and 0.56 grams Ag/cc compared with 
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Figure 24 


Plot of Ln D versus 1/T for the Incremental Thin Film 


Couples Annealed for 24 Hours. 





shit 


the numerical solution with the corresponding times entered 
into equation (19). 

Grain size determinations, using the transmission mode 
of the electron microscope, were made on each couple... ~The 
grain size varied between 1,000 and a OIGK0 K on all couples 
observed. No variation across each couple, including the 
diffusion zone was observed. Figure 25 shows photomicro- 
graphs of the grain distribution in the 200 and 350°C 
couples. Figures 26 through 28 are photographs of the 
surfaces of selected incremental couples Showing the change 


in surface structure with increasing deposition temperature. 


Ag-Au Couples 


Table 5 is a summary of the experimental test conditions 
of Runs 201 - 223. Each run involved the deposition of a 
pure gold film layered side by side with a silver film of 
equal thickness. 

Runs 202 through 212 have comments in Table 5, listing 
reasons for the author's inability to use them for quantita- 
tive results. The problem occurring in the initial runs 
(202 - 206) was film breakup (dewetting) of the gold side of 
Eiewcouple. Along with this, there was some interesting 
Pareicle Shaping of the resultant discrete gold islands 
observed, such as smooth equilateral triangles, like those 
observed [40] in epitaxial gold films on [111] surfaces of 


M,S,, and particle shapes observed by Sundquist [41]. 





a 
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Pigeine 25 


Transmission Electron Micrograph of Section of the Surface 
Oreinsupportedeliany Film, 23,000~x 
A. Gold Side of Incremental Thin Film Couple 
Annealed at 200°C for 24 Hours. 
B. 80%Au-20¢Ag Side of Incremental Thin Film 


Couple Annealed at 350°C for 24 Hours. 
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BO Au - 2OAg 
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Eoure a2 6 
Scanning Electron Photograph of the Surface of the 
Unannealed Incremental Thin Film. Couple (Run 227) Showing 
Interface Region Film at 45° Tilt to the Electron Beam 


10,000x. 
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Scanning Electron Photograph of the Surface of the 
Incremental Thin Film Couple Annealed at 250°C for 24 Hours, 
Showing Interface Region. Film at 45° Tilt to the 


Electron Beam. 10,000~x. 





Se) 





Risure 26 


scanning Electron Photograph of the Surface of the 
Incremental Thin Film Couple Annealed at 416°C for 24 Hours, 
Showing Interface Region. Film at 45° Tilt to the 


Electron Beam. 10,000~x. 





wip ic 5. 
Temperature 
oils 
Anneal* 

Run ioe) 
201 aS 
202 675 
203 SBS) 
204 500 
205 675 
206 500 
207 428 
208 416 
209 416 
ZrO 416 
Z11 416 
Zi 2 416 
JEM 416 
Zia 416 
ZS 416 
ZO 416 
Zi, 416 
218 416 
219 416 
27) 300 
aia Al 300 
CeZ 300 
275 300 
a ag OD 

geo OA 


+ 
i 


not measured on Dektac 


6 


summary of Experimental Conditions of Runs 201 


three 223. 


Time 
of 


Anneal 


(hrs) 
0.0 


Fo-..00 


24. 
48. 
48. 
48. 
48. 
48. 
90%. 


eee Go. > C5 SO Se Se 


Vice 
24. 
48. 
2 


48. 
ZA. 
24. 
ance 


dO 
ee eG 2 a 


Film 


ihe kies s 


(A) + 
1000* 


1000# 


1000+ 
1000* 


1080 


1000 
1000# 
1550 
1500# 
15004 
15004 
15007 
1000 
1000 


1000 
1520 
1540 
a0 


1490 
1560 
1500# 
1470 
1490 


(Ag-Au) Couples 


Comments 


used as time = 0 profile 
for Ag-Au couples. 


excessive breakup 
(dewetting) of gold 
Pportionsemcouple 


Same as 202 


Partialedewetting of 
gold side; no interface 
COntacE 


excessivemcurtace rough- 
ness of gold and silver; 
partial dewetting 


Same as 205 

no interface contact 
wide region of mixing 
no interface contact 
same as 209 

Same as 209 

same as 209 

run 201 annealed 


run 201 annealed; 
dewetting of gold 


run 201 annealed 


no interface contact 
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Three (3) runs, 213, 214 and 215 were made, using 
sections of the 1,000 A unannealed couple (Run 201), at 
megeceior times of 48, 72 and 96 hours. The formation of 
a porous zone, running parallel to, and on the silver side 
of the interface, was observed in the couples annealed for 
48 and 72 hours. In addition, a zone of protrusions was 
mmeerved On the gold side of the interface. The porous 
zone was so extensive that complete film separation occurred. 
The width of the porous zone was slightly larger on the 72 
hour couple. In the couple annealed for 96 hours, the film 
separation due to the porous zone was evident. However, 
dewetting was observed on the gold side of the film, with 
particle shaping like that described above for films 
deposited at higher temperature. The silver side of the 
couple also exhibited extensive pitting, and investigation 
on the scanning electron microscope revealed a high degree 
of surface roughness. Subsequently, no quantitative study 
G@eedattusion was undertaken for these runs. 

Runs 216 through 219 involved the deposition of 1,500 A 
thin film couples at 416°C, for annealing times of 2, 12, 

24 and 48 hours. Again, the formation of the porous zone 

and surface protrusions were observed on the silver and 

gold side of the original interface, respectively. Figures 
29 through 32 show photomicrographs of the respective inter- 
face region of the four (4) couples. Figure 33 1s a scanning 
electron photograph of the interface region showing the 


porous zone, and the associated surface protrusions. 





98 


Av 





“| top. be 


Pacures.g9 
Optical Microscope Photograph of Interface Region of Ag-Au 


Thin Film Couple Annealed at 416°C for 2 Hours. 1080~x. 
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Paieure 250) 
Optical Microscope Photograph of Interface Region of Ag-Au 


Thin Film Couple Annealed at 416°C for 12 Hours. 1080x. 
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Paleo eon 
Optical Microscope Photograph of Interface Region of Ag-Au 


Thin Film Couple Annealed at 416°C for 24 Hours. 1080x. 
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Figure 32 
Optical Microscope Photograph of Interface Region of Ag-Au 


Thin Film Couple Annealed at 416°C for 48 Hours. 1080~. 
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Scanning Electron Photograph of Interface Region of Ag-Au 
Thin Film Couple Annealed at 416°C for 24 Hours. 


Reimeaeee 5 lilt) to the Electron Beam. 3,000~x. 
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An investigation of the chemical composition of the 
surfuce protrusions was undertaken by measuring characteris- 
tic x-ray lines generated from the electron beam of the 
Beanning clectron microprobe. For a counting time of 200 
seconds, the counts of the La lines of the silver and gold 
remained constant with the electron beam, on and off, 4 
different surface protrusions. The silicon counts decreased 
to one-half their original value when the electron beam was 
On a surface protrusion. 

Characteristic x-ray measurements were also taken for 
Bums 216 - Z19 as a function of lateral distance on the 
Seectron microprobe. <A computer program similar to that 
used for the incremental couples, converted the intensity 
Gomnts Co units of mass concentration. The direct computer 
Output of the couples are shown in Figures 34-37. 

Because of the formation of the porous zone and associ- 
ated protrusions, there was film thickness variation as a 
mumetton OL lateral distance, thus the concentration-distance 
profiles were suspect. 

A scanning electron microprobe investigation of the 
couples, however, indicated the presence of silver on the 
gold side of the porous zone extending away from the inter- 
face approximately the same distance as indicated by the 
concentration-distance profiles. No counts of gold were 


observed on the silver side of the porous zone. 
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The dip in the concentration value observed in the 
concentration-distance profiles was taken as a representa- 
tive of the center of the porous zone. A computation of the 
area under the concentration profile from the center of the 
dip on into the gold region (increasing relative distance) 
was made on each profile. The area was determined using a 
computer program written to numerically integrate, employing 
Simpson's rule. The results are summarized in Table 6. 

A computer program similar to that of Appendix E was 
Piecen to treat the peak in concentration observed in the 
4 couples. Run 219 (Figure 34), the couple annealed for 2 
hours, was treated as the initial profile (time = 0) and 
the values of concentration on the gold side of the porous 
meme was fed into the computer program along with a compli- 
mentary set of image concentration values to make the initial 
profile an even function. The reduction in peak height of 
the 12, 24 and 48 hour annealed couples was compared with 
the numerical solution for times of 10, 22 and 46 hours 
GFespectively. the numerical solution to the concentration 
“nulse" of Run 219 is drawn on each concentration profile 
AOweaneappropriate value of DBD. D values of 1, 5 and 9x10 *cm?/ 
sec fit the experimental profiles for the 12, 24 and 48 hour 
ammeal, respectively (Figures 35-37). 

Runs 220 - 223 were 1,500 A thick thin film couples of 
Ag and Au annealed at 300°C for 12, 24 and 48 hours. The 


formation of a porous zone and associated surface protrusions 





Ps 


Table 6. Comparison of Area under Concentration Curves from 
the Edge of Porous Zone into the Pure Gold Region 


of the Couple. For Runs 216-219, Annealed at 416°C 


Time of Area 
Microprobe Anneal Ceonc, Average 

Run Scan No. (hours) mveT ons) Area 
219 1 2 61.9 58.4 
219 2 2 54.8 
218 1 i 62a 60.2 
218 Z he 58.4 
216 1 24 54:20 oe 
NS E 24 49.3 
ZAG i 48 a0 Sd ed 


a, 2 48 S305 





TEs: 


were again evident. The width of the porous zone of the 12 
hour annealed couple was similar to that of the 2 hour 
annealed couple at 416°C. The porous zone enlarged with 
time, and film separation was evident in the couple annealed 
for 48 hours (Run 220). X-ray measurements as a function 

of lateral distance were made on the electron microprobe 

and the data were converted to units of mass concentration 
as described above. However, the porous zone of the 12 hour 
annealed sample was so small in extent that the dip in con- 
centration was not evident. Therefore, the semi-quantita- 
Mme technique to determine the rate of silver penetration, 


as described above, was not attempted. 








aoe? 


DISCUSSION 


Thin Film Standards 

The linearity (high degree of correlation) observed 
between the microprobe intensity and film thickness for all 
the data exhibited in Tables 1 and 2 suggests that the thin 
film approximation, intensity is proportional to mass thick- 
ness, 1s a valid assumption throughout the entire composition 
range of the gold-silver system, for the film thickness 
studied. The reader is cautioned, however, that the numeri- 
cal measurements are only valid for gold and silver films on 
a silicon substrate with a 1,500 RK oxidized layer on the 
Suri ace . 

The statistical equations used in the analysis of thin 
film standards, and subsequent work on thin film couples, 
are shown in Appendix F. 

The maximum deviation from linearity occurred with the 
2 UUW A films as should be expecucea,.cspecially in the casc 
of pure gold films where I/V, ik was approximately 0.6. The 
ani cehative Crror, #2o,01n the cmpirical conversion 
parameter, was observed in the measurements taken off the 
Silver spectrometer for films between 500 and 2,000 A. 
Peealsemorethe small relative error, One conversion param- 


eter was used. The parameter was determined from a linear 
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reast squares fit of all the data, representing the relation- 
ship between measured relative intensity and weight fraction 
composition for films of thickness between 500 and 2,000 Ne 
Again, the largest deviation from the straight line fit 
Becurred for the 2,000 A films of high gold composition, 
where the thin film approximation became suspect. 

The relation between an error in the empirical conver- 
Slon parameter and errors in the calculated composition value 


1S given by the usual expression 


Peace 
and from equation (2) 
al Geen : 
a CUl=C) 7/4. (25) 
mies error in composition is 
AC = U(1-C) Aasvea (24) 


) 


Bd 


The maximum absolute error in composition is one-fourth ( 
Stecie relative error in the conversion parameter, or 
ie25)(.25) Which is approximately 6 weight 4%. 

This, however, is based on the assumption that the com- 
position of the standards is known to within negligible 
limits of the error of the calculated value. 

The determination of the composition of the thin film 
standards used in the experimental work of this author was 


based on the assumption that the thin film density in both the 


pure and alloy states was that of the reported bulk value 





Vee) 


discussed in the Background and Literature Survey section. 

Hartmann [{42] obtained bulk values for Au and Ag at 
thicknesses down to 200 A. However, Blois and Reiser [43] 
found the density of Ag films to approach bulk values at a 
thickness of 1,200 A. The data of [4 sivare suspeces. / |; 
for silver, which readily forms a surface oxide layer, 
mowering the effective density of the film. 

This author, in the initial stages of determining the 
Orientation relationship between the mass deposited on the 
monitoring crystals and that of the substrate, made three 
determinations each between the frequency shift of the 
crystal and measured thickness (Dektac) of pure gold and 


emiver films. Using the relationship, 
Af = Bpt , (25) 


where Af = frequency shift of the crystal 


t = film thickness 


SUboerate sens ual solmcerOnlentation parameter , 


B 


0 Heol Cl © il Sedey, 

a linear relationship was observed between Af and the 
measured value of film thickness. Solving equation (25) 

OT Pag and Pry respectively, a bulk value was obtained 
within the relative error of the measured film thickness. 

If the standards are assumed to be of bulk density, then the 


major source of error is the measurement of film thickness, 


+ 50 A. The error in composition can be written 
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Mee t 
AC(atomic 3) = ——S& Ag _ (26) 
At, At 
a Au 
CAg CA 
which has a maximum value of + 10 atomic % for film thickness 


1] 
of 500 A and at each extreme of the composition range. 


Ag-Au Couples 


The experimental observations of film breakup on the 
gold side of the thin film couple suggest that the mechanism 
of atomic reorientation manifests itself as a function of 
time and temperature, and can be attributed to an annealing 
effect. The agglomeration of the gold side of the couple 
observed in Runs 202 - 206, at substrate temperatures above 
500°C, was only observed at the conclusion of the anneal. 
The initial condition of the film, immediately after film 
deposition was in doubt and the observed effects, when com- 
pared to that of Run 201, could only be attributed to the 
substrate temperature. However, in the series of runs at 
milo -€ (Runs 213 - 215) the 1,000 A film deposited at room 
temperature (Run 201) was used, and the following anneal led 
to agglomeration of the gold after 96 hours (Run 214). 

If a polycrystalline material is annealed, grooves 
develop at the intersection of the grain boundaries with the 
metal-vapor surface. If it is assumed that the grain bound- 
aries do not move during annealing, then it can be expected 
that the groove depth (d) deepens with time. Mullins' [44] 


well known analysis of grain boundary grooving kinetics was 
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applied to the thin films of the present work, thus 
(0.486)y,, nts | 
> “al oan ie) 
where d= depth of groove 
= cpcCliilceinee energy of the grain boundary 


cus Specific free energy of the surface 


De = surface diffusion coefficient 
k = Boltzmann's constant 

T = temperature of anneal (°K) 

t = time of anneal 


N = number of surface atoms/unit area 


i) 


volume per atom 

In order to relate the experimental observations of the 
gold film breakup, a two stage model was proposed. In the 
initial stage there would be groove deepening to the sub- 
menace, the depth of the groove equal to the original film 
thickness. Secondly, the hole grows, with critical size 
agglomerates (adjacent to the holes) growing at the expense 
of the surrounding metal. 

If equation (27) is evaluated at 416°C, (temperature of 
anneal of couple where gold agglomeration occurred after 96 


hours) using the following values, 


vent 1450 from [49] 
ane 450 from [49] 
Q = 1.7 x 10°77? cm? 


eee eS wena 


Foe 
i 
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te = > 0s” cm /Seemhor Auvat 416°C from [51] 


ee Oe 07S Sei ercemior Neat 416°C from [51] 


The depth of the groove would equal 1,000 A (film thickness) 
for annealing times of 62 and 2.5 x 10* hours for gold and 
silver respectively. 

In the case of gold, we would expect groove deepening 
to the substrate after 62 hours. Following this, the anneal 
enables the film substrate system to lower its free energy 
by lowering its interfacial free energy. The motive force 
for the growth of an existing hole by surface diffusion is 
the surface tension vector sum of the film-substrate, sub- 
strate-vapor and film-vapor interfaces. The flow of migrating 
mpems 15 due to the difference in chemical potential at regions 
of differing curvature. 

Presland et al., [45 and 46] in a series of two articles 
applied the theories of Mullins [44] to his experimen ccal 
observations in annealed thin films of silver. He showed 
that following an induction period (hole deepening to the 
Substrate with subsequent hillock growth) that further hillock 
growth followed, with critical size agglomerates growing at 
the expense of surrounding areas. The continuing hole growth 
resulted in a decrease in metal-substrate interfacial area, 
resulting in the separation of agglomerates (film breakup). 

Brandon and Bradshaw [47] showed that the surface self 
Giriicion coertricient Of Silver iS increased in the presence 


of oxygen. Presland could only produce agglomeration in his 
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high vacuum deposited silver films when the silver vapor- 
metal interface was exposed to oxygen pressures greater than 
Ms torr. Below Po, = 10-*’torr, he reported groove deepen- 
ing with subsequent hillock growth. However, the catchment 
area for the critical size agglomerates increased rapidly 
with decreasing Po,? Suggesting changes in surface tension 
and the surface self diffusion coefficient with oxygen 
partial pressure. In comparison, gold agglomeration has been 
Observed [48] in unsupported thin films in high vacuum. 

The equilibrium shape of a particle on a substrate is 
determined by the vector sum of the interfacial free energies 
of the film, vapor, substrate system. If the observations of 
[46 and 47] are assumed to be operative for the gold film, 
then in conjunction with the sufficient groove deepening, as 
Predicted by equation (27), the vector sum of the interfacial 
meee Energies iS in accord to predict hole growth and subse- 
qaent agglomerates. 

In the case of silver, the reduction of the film-vapor 
interfacial free energy, due to the presence of oxygen on 
the free surface of the silver as suggested by Presland [45, 
46], would not be sufficient cause for hole growth. 

More likely, the increase in the surface self diffusion 
coefficient due to oxygen pressure would reduce the time 
necessary for groove deepening to the substrate. Following 
groove deepening, a reduction in the substrate-vapor inter- 
facial free energy due to oxygen pressure would be favorable 


for hole growth followed by film agglomeration. 
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The observed faceting on the gold film agglomerates 
weuld be attributed to the reported [41] surface anisotropy 
of gold. Gjostein [50], who provided some of the impetus for 
the studies of [45 and 46] in an earlier paper [51], is 
currently developing a manuscript to possibly improve on the 
studies of [45, 46 and 47] and should provide interesting 
reading when released. 

To avoid the agglomeration of the gold film, further 
runs (216-223) were made at temperatures below 500°C and for 
annealing times only up to 48 hours. In each of these runs 
the porosity observations of bulk investigators with pure 
Silver and gold sandwich couples were similarly reproduced 
Baume thin film couples. The porosity increased with time 
and temperature of anneal. However, in contrast to the bulk 
Counterpart, physical separation of the films occurred in 
the lateral direction. 

The physical separation of the two regions of the film 
leads to thickness variations of the couple around the porous 
zone. The concentration profiles of Figures 34 to 37 were 
determined by assuming a uniform thickness across the diffu- 
Som. ZONE . Te neeones in the region of the porous zone, the 
concentration profile is in error. 

Seeunmdily, the physical Separation of the couple severely 
limits the atomic mobility of the two species. When this 
Mectrs. the quantity of pold and silver which had previously 


diffused into the silver and gold respectively, remains fixed 
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with time. The observation of the similarity of the areas 
imeer the concentration profiles of the 2, 12, 24 and 48 hour 
anneals confirms this conclusion. The smaller area observed 
mimeme 48 hour anneal (Run 217) can be attributed to the 
inability of the electron microprobe to detect the small 
concentration of Ag above background at distances far from 
the porous zone, and the corresponding increase in error of 
the empirical conversion parameter with low values of solute 
momcentration. 

The mass of silver remaining fixed with respect to time, 
eeter the formation of the porous zone, creates a set of 
boundary conditions different from that of an initial t = 0 
profile. The series of photographs of the four anneals shown 
iMmimerreures 29-32 indicate that physical separation of the 
couple is almost complete after the 2 hour anneal (Figure 29). 
On cither side of the porous zone the film thickness is 
uniform, therefore the concentration profile of the 2 hour 
famed! can be taken to represent the new initial profile 
ee = 0) and the 12, 24 and 48 hour anneals to represent 
weimes Of t' = 10, 22 and 46 hours respectively. 

The thin strip approximation was applied to the concen- 
tration "pulse" of silver observed on the gold side of 
Bomolcezone On the t” = Oh protile. the center of the “pulse 
peak" indicated on Figure 34 was taken to represent the edge 


of the porous zone shown in Figure 29. 





126 


ine exact solutvonmpe2|eto the initial (t* = 0) concen- 
tration profile would be to represent the profile by f(x) 
and solve using a heat flow analogy. The porous zone would 
represent an impenetrable barrier (insulating medium). 
A numerical approximation to the exact solution was used, 
Incorporating the thin strip approximation. The concentra- 
tion "pulse" was made into an even function by reading into 
a computer program the image concentration values as shown 
in Figure 34. The thin strip approximation predicts a 
reduction in peak height with points of constant concentra- 
tion moving laterally parabolically with time. The resultant 
fit to the t' = 10, 22 and 46 hour profiles shown in Figures 
feeto 57 lends validity to the above approximation. 

ine =peductiongim ) as the time of anneal increases can 
again be attributed to the errors in detecting small solute 
(Ag) concentration values and the rising error in the con- 
version parameter with decreasing solute concentration. 
Also, the errors in the thin strip approximation decrease 
as YDt increases relative to the assumed (0.8 microns) strip 
thickness (lateral direction). The value of D therefore, 
is probably somewhere between the results for the t' = 10 
and 22 hour anneals. Averaging the two results gives a 
D of 7 x 10-11 cm?/sec. which closely represents what was 
reported for the incremental couples at 416°C. 

An interesting observation was that of the secondary 


zone Of porosity observed in the photographs for the 24 and 
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48 hour anneals (Figures 31 and 32). This could possibly be 
attributed to the new set of boundary conditions operative 
after the 2 hour anneal. The mass flow of silver in the t'! 
framework could be viewed as a thin film of gold layered 
side by side of an alloy film of say 50 atomic $ silver with 
only a 30 micron extent (5 width of concentration "pulse"') 
in the lateral direction. The new set of boundary conditions 
would also shift the position of the original interface and 
the corresponding position of the maximum vacancy flux 
Pragient [53]. 

This reasoning implies that the vacancy flux gradient 
must be a strong driving force for the nucleation of voids 
by an assimulation of excess vacancies. On the basis of 
energy considerations it would appear logical for the 
vacancies accumulating in the subsequent mass depletion zone 
to continue to annihilate themselves on the already formed 
voids in the original porous zone. 

ihe X-ray investigation on the scanning electron micro- 
scope of the surface morphology on the gold side of the diffu- 
Sion zone confirmed the suspicion that the morphology repre- 
sented positions of increased film thickness or surface pro- 
GEusions. 

Lomwitheeche PpOrostty 1Ormatlon, nNowin-depth investigation 
of the formation mechanism of these surface protrusions, or 


if the reader will allow, surface hillocks, was undertaken, 


for it was not the purpose of this work. 
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llowever, some speculations will be made from the previ- 
ous work of bulk investigators in the Ag-Au system and those 
studying thin film annealing phenomenon. 

several investigators, [54, 55 and 56], studying thin 
film annealing phenomenon have reported hillock growth in 
thin films, including Au, and attributed it to a mechanism 
mumestress reduction. The driving force in their studies 
was the differential stress between the film and substrate, 
developed by the AT between the deposition and annealing 
memperature. In the present work, deposition at the anneal- 
mie temperature (excepting Runs 212 - 215) eliminated this 
as a driving force. However, Barnes [30] and Balluffi [32] 
reported strain in the diffusion zone of bulk sandwich 
Geuples Of the Cu-Ni and Ag-Au systems respectively. fhe 
memaction of strain in the diffusion zone and/or dislocation 
climb creating planes of vacancies similar to that proposed 
Dy Ruth [35] in the silver-gold system, could account for 


the observed hillock growth. 


Incremental Couples (80 Au-Au) 
As in the bulk studies of Johnson [38] with sandwich 
Souples with small initial concentration differences, no 
porosity or hillock formation was observed in the diffusion 


zone of thin film couples with an initial 20 atomic % con- 


centration difference. 
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The transmission electron microscope studies of the 
mouples’ grain size led to the conclusion that the grain size 
mev00 - 2,000 A) of the films was independent of the anneal- 
ing temperature, and film composition over the working range 
of the couples. The grain size of Ag films deposited at 
room temperature has been reported to be as low as 50 A. 
However, with the high deposition temperature of the present 
Study, the grain size is expected to be larger, due to the 
increasing surface mobility of adatoms and clusters during 
deposition. 

The surface roughness, which was uniform over each 
individual couple, increased with annealing temperature as 
feieeonoOwn in Figures 26 to 28. As the surface mobility of 
adatoms increases (higher substrate temperatures), the con- 
demsation can occur preferentially at the surface concavites, 
and thus tend to smooth the surface. However, as was 
reported [41] with gold and experimentally observed in the 
present study, anisotropy of the surface energy will cause 
the condensate to grow preferentially along certain crystal 
faces, enhancing the resultant surface roughness. This is 
strongly evident in the couple annealed at 416°C shown in 
Repure 28. 

In each of the concentration profiles, the variance in 
the slope was within .005. The value of .005 was within the 
statistical error of each slope. The error was determined 
from the experimentally observed variation in slope between 


successive scans of the microprobe across the diffusion zone 
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teoiaiicrent portions of the couple. The three values of D 
demerminced at 416°C for annealing times of 12, 24 and 48 
mours lent credence to the use of the thin film approximation 
which predicts parabolic movement of atomic planes of 
constant composition with time. 

The values of D determined for the 24 hour anneals at 
memperatures of 200, 250, 300, 350 and 416°C exhibited 
linearity on a In D versus 1/T plot with correlation 0.99. 
This indicated a single atomic transport mechanism to be 
dominant over the working temperature range of the couples. 
The value of activation energy was almost ¢ of that reported 
romepulk diffusion in Ag-Au couples. 

Mimic rnin ti imvecouple — three, possible atomic transport 
mechanisms are considered to be operative: Transport of 
matter along the grain boundaries and the substrate-film 
interface, through the crystal lattice, and along the f1ilm- 


vapor interface. The apparent diffusivity can be written as 


D = D S (28) 


aoe cp & + D, (1-g-s) + D 


S 


where Dep oral DoOuUnatnyedinr US TOMmecoctfi client 


lattice. diffusion coerLicient 


oO 
Ul 


L 
De = Shr tace Gdirtuls1omecoet1lerent 
g = fraction of time a diffusing atom spends on a 
grain boundary or along the substrate-film 
interface 
S = fraction of time @ diftusing atom spends on the 


surface. 
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LOGetncwGase Of aupanery dlloy, all atoms are net equiv- 
pent SO the total time a diffusing atom spends on a grain 
boundary or surface is equal to fractional number of sites 
which are in or on each respectively, times its atomic 
meaction. 

To estimate g, consider a sample made up of "'spheres" 
that touch each other completely. The total number of atoms 
in a "sphere' is 4/3mr°/2, where 2 equals the atomic volume. 
The total number of atoms in the grain boundary equals 


eeeya!: Therefore, 





g = (atm. grac.) [A> [atts] = (atm. frac.) (32) 


fOr gifs = 2.2 A 
r = 1,000 A 


Oil 


atm. frac. 


Therefore g = 10°? 

To estimate s, consider a surface monolayer of thickness 
2 A. For a film thickness of 1,500 A and atomic fraction 
ies = 10°. 


We can rewrite equation (28) as 


= = 5! ~ & 
Dapp 7 20 Dey 2 OO (29) 


The activation energy for chemical diffusion through the 
lattice is approximately 40 kcal/mole, which is not in accord 
with the activation energy reported in this study. Over the 


temperature range of this present work the contribution to 
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Eire apparent diffusivity using the bulk value of diffusivity 
reported by Johnson [38] is negligible. 

If surface diffusion is assumed to be dominant, then a 
Similar discrepancy exists between the reported activation 
energies for surface diffusion and those of the present work. 
If we consider the terrace ledge kink model of a solid-vapor 
interface developed by Kossel [59], then the activated 
process for transport along a (111) FCC surface involves two 
SIECDS . 

The activation energy is considered to be the sum of two 
contributions, a change in enthalpy necessary for formation 
of the diffusing species and secondly, the change in enthalpy 
necessary for motion. The adatom in general carries the 
surface flux. It is the most weakly bonded atom to the 
pierace, and on a FCC surface is characterized by 3 nearest 
neighbor bonds. An adatom is formed by an atom escaping from 
a kink site in a monatomic ledge (6 nearest neighbor bonds) 
to a site at the ledge (5 nearest neighbor bonds) and then 
to an adatom site. This involves the breaking of 3 bonds. 
The motion of an adatom involves the breaking of one bond. 
The activation energy therefore, can be looked upon as a 
process involving the breaking of 4 nearest neighbor bonds. 
The heat of sublimation of a FCC metal is the enthalpy change 
when a mole of atoms each atom breaking 6 bonds, escapes to 
BRHOMMapOn State. Therefore, we can estimate the activation 


energy per mole of atoms as 2/3 AH For gold and 


Sub lama tion 
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Silver, the heat of sublimation is approximately 80 and 60 
kcal/mole respectively. This high value of activation 
energy 1S in accord with what has been reported for gold and 


Silver where T/T 1s greater than 0.77. 


Ment. pit. 

It has been suggested by Gjostein [51] that for studies 
at low temperatures (T/Tm < 0.77) the activation energy is 
approximately one-fourth the heat of sublimation. One 
study [60] of grain boundary grooving in gold reported an 
activation energy of 17 kcal/mole for T/Tm = 0.57. 

bso rder to be an accord with the presently accepted 
ferrace ledge kink model, it must be postulated that the 
ferracce 1S cOvered with a large fraction of growth sites 
available for surface transport. The activation energy can 
then be assumed to involve the breaking of one bond. This 
model is for the most part unsubstantiated and sufficient 
experimental work has not been done in the present study to 
Support Or Challenge this model. 

If we consider grain boundary diffusion to be the domi- 
nant mechanism, then we again are faced with interpreting the 
low value of activation energy experimentally observed (10,900 
cal/mole). Hoffman and Turnbull [36] reported an activation 
of 20,200 cal/mole for a study of grain boundary diffusion of 
radioactive silver in a silver bi-crystal. If we just take 
into account the grain size of the film, the ratio of 
/Dop equals 0.01, 0.15 and 0.40 for temperatures of 200, 


D 
app 
300 and 416°C respectively. Let us assume the results of 
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Hoffman and Turnbull [36] as representative of a lower limit 
on the degree of disorder in a grain boundary. An upper 
limit would be to consider the grain boundary as in the 
liquid state. Activation energies for liquid self diffusion 
in silver have been reported [58] between 7 and 8 kcal/mole. 
The ratios of the apparent diffusivity to diffusion down 
"liquid" (highly disordered) grain boundaries are greater 
than 1 for all annealing temperatures. If we assumed that 
the grain boundaries between the grains in the thin film 
couples are highly disordered as the photographs of the film 
mrainm structure in Figure 25 and Figures 26 through 28 of 
the film surface suggest, then we would expect that the 
activated energy for transport along the grain boundaries 


to be between the upper and lower limit discussed above. 
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CONCLUSIONS AND RECOMMENDATIONS 


Conclusions 

The study of lateral diffusion in thin film couples of 
gold and silver can be investigated using one of the non- 
Steady state solutions to Fick's second law. The thin strip 
approximation, which in the limit of strip thickness = 0, 
is more commonly applied to a pair of semi-infinite solids 
faledlisoO be applied to a thin film of constant thickness 
where the initial profile is experimentally determined. 

The use of Zeibold's empirical conversion parameter, 
more commonly applied to the determination of composition 
in the bulk, is valid for thin films of the gold-silver 
system up to 2000 A in thickness. Up to 2000 A the thin 
film approximation that intensity 1s proportional to mass 
thickness is valid in the gold-silver system. 

Studies of lateral diffusion in Ag-Au couples exhibit 
surface morphology changes similar to that reported in bulk 
couples. The mechanism, operative in the thin film diffusion 
couple giving rise to the formation of a zone of porosity and 
protrusions on the silver and gold sides of the interface of 
Mremeouple respectively, can be attributed to the conclusions 
reported by bulk investigators. The maintenance of continuity 


of the silver and breakup of the gold sides of the thin film 
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couple follows the predictions of Mullin's theory of grain 
boundary grooving. 

The use of incremental couples with small initial con- 
centration differences eliminates the formation of porosity 
and hillocks as experimentally observed in the bulk studies 
of Johnson. The maintenance of continuity in film thickness 
across the diffusion zone allows the determination of diffu- 
Sivity as a function of temperature. Using the standard 


fecmecquation for an activated process the apparent diffusiv- 


ity was written as D = 1.9 x 107’ exp [12,00 cm*sec~? where 
D, = 1.9 x 10°’ and the activation energy is 10.9 kcal/mole. 


In addition a value of diffusivity was determined on the 
Ag-Au couple which had film separation by considering the 
porous zone as creating a new set of boundary conditions for 
Bie wexpcrimentally determined initial profile. The value of 
Mperusivity was in accord with that observed in the incre- 
mental couples. 

Pecomparisonm Of Che possible Cransport mechanisms 
Operative in the thin film couples between 200 and 416°C 
Suggests that grain boundary diffusion is the dominant 
mechanism. In order to resolve the low value of activation 
energy reported in both the incremental couple and the Ag-Au 
couples an assumption that the grain boundaries were highly 


disordered must be made. 
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Recommendations 

[In order to make a stronger case for the mechanism of 
grain boundary diffusion, incremental couples should be made 
with varying rates of deposition in order to vary the 
mesultant grain size. larger grains with a rising value of 
activation energy in accordance with the previously reported 
bulk values of grain boundary diffusion would lend more 
credence to the disorder assumption made in this study. 

A study of incremental couples with larger concentration 
differences could prove interesting. The magnitude of the 
initial concentration difference necessary to cause the 
formation of porosity could produce some interesting con- 
clusions as to the operative driving force for the formation 
@eethe porosity and hillocks. 

Because of the small amounts of material present in 
thin films and short annealing times, a study of lateral 
fer tusion using the scanning electron microprobe could prove 
valuable. The small beam size would allow greater resolu- 
tion, in the lateral direction, and in the relationships of 
surface morphology to chemical composition, than in the 
electron microprobe. 

Merollowin study (on ene sapplication of the empirical 
GOnversion parameter of Zeibold to other thin film systems 
could prove to be a valuable tool in film composition deter- 


minations. 
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The study of non-miscible systems using the experimental 
deposition technique of the present work should prove 
superior to the overlap technique used in previous work. 

The elimination of the overlap region removes the necessity 
for diffusion normal to the substrate to homogenize the 
Overlap region and more closely approximates the initial 
conditions of the bulk sandwich couple. The method of 
deposition used in the formation of the incremental couples 
could be applied to the study of diffusion in non-miscible 


Systems in single phase regions. 
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APPENDIX A 
A STUDY OF LATERAL DIFFUSION IN 
THIN FILM COUPLES OF BISMUTH AND TELLURIUM 


This author's initial investigation of lateral thin 
film diffusion took place in the bismuth-tellurium system. The 
separation between the razor blade and the substrate was 
achieved using a differential screw driving mechanism [9]. 

The adjustment was done manually and resulted in regions of 
mixing on the order of 10~° meters. 

All Bi-Te thin film couples were 1,000 A thick and 
deposited at room temperature. The anneals were carried out y! 
in two separate environments. Vycor encapsulated tubes at a! 
residual gas pressure of approximately 50 microns were used 
for anneals at 1, 3, 10, 30 and 96 hours. In addition, 
anneals in high vacuum at pressures less than 10°’torr were 
carried out for times of 1, 3, 10, and 30 hours. 

The initial profile was determined on the electron 
microprobe, on an unannealed thin film couple. The region 
of mixing was 84 microns. From optical observation, and 
the electron microprobe intensity data, 71.5 and 12.5 
microns were attributed to the bismuth and tellurium overlap 
region respectively. The interface (position of razor blade) 


separated the two overlap regions. 
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In the series of anneals in the vycor tubes, the follow- 
ing observations were made: A peak in tellurium intensity 
was observed to move through the region of mixing towards 
the pure bismuth side of the couple into the tellurium over- 
lap region. The anneal at 10 hours signaled the completion 
of the shift of the tellurium "pulse," which now was entirely 
in the original tellurium overlap region. In the bismuth 
wer lap region (tellurium side of original interface) a 
reduction of tellurium on the pure tellurium side of the 
couple was observed, followed with an increase in tellurium 
myer the entire bismuth side of the couple. The pure 
tellurium film thickness originally at 1,000 A was reduced 
to 200 A. 

The observations of the couples annealed in high vacuum 
were the same for 1, 3 and 10 hours. However, the 30 hours 
anneal exhibited no significant decrease in tellurium 
@mekmess on the pure tellurium side of the couple or 
subsequent increase in tellurium intensity over the bismuth 
side of the couple. 

The bismuth intensity profile remained immobile for all 
annealing times in both environments. X-ray fluorescent 
aiabysis Using Cu Ka radiation Confirmed the presence of the 
intermetallic compound Bi2zTe3 in the region of mixing for 
annealing times greater than or equal to 3 hours. After 10 


hours no further growth of the intermetallic compound was 


observed. 
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The mean free path of ambient gas atoms at 50 microns 
m0. 5 cm. Tye fate Of ettveemevaporation of tellurium and 
ememuth at 200°C is 0.2 and 2 x 10-° monolayers/sec.., 
respectively. Therefore, mass transport by the vapor phase 
is a probable mechanism in the couples encapsulated in the 
vycor tubing. This would explain the increase and decrease 
in counts of tellurium over the bismuth and tellurium sides 
of the couple for the 30 and 96 hour anneals in the vycor 
Capsules. 

The formation of the intermetallic compound occurs in 
both environments, ruling out a vapor transport mechanism, 
especially unlikely in high vacuum at 10-’torr where the 
mean free path is 450 meters. 

The relative immobility of the bismuth profile, and the 
Meving tellurium "pulse" for the anneals at 1, 3, and 10 
hours in both environments led the author to conclude that 
the tellurium was the only active (within the time frame of 
Bae author's study) species in the mass transport process. 

Let us consider the region of mixing to be made up of a 
Series of thin strips side by side in the ea direction. 
Because of the simultaneous deposition process, each strip 
will be homogeneous in composition in the direction normal to 
the substrate. The variation in strip composition in the 
lateral direction should be from pure tellurium to pure 
bismuth across the region of mixing starting from the pure 


tellurium side of the region. Therefore, one strip in 
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Original region of mixing should be of the stoichiometric 
composition of the intermetallic compound Bi2Te3. If we 
arbitrarily assign the mass ratio of bismuth to tellurium 
(m,,;/mr.) in this strip of stoichiometric BizTe3 equal to l, 
Mien Strips On either side would have a mass ratio < 1 and 

> 1 towards the pure tellurium and bismuth side respectively. 
Mass transport of tellurium in the direction of the pure 
bismuth side of the couple would increase and decrease the 
respective mass ratios towards 1, on the tellurium and 
bismuth sides of the initially stoichiometric strip. The 
region of mixing would now be a series of strips of Bi2Te3. 

However, after the 10 hour anneal which signals the 
conclusions of the mass transport process, the region of 
mixing, and now zone of intermetallic compound, appears to 
remain stationary for further times. If the compound forma- 
tion was diffusion controlled, a parabolic growth of this 
zone with respect to time, would be expected. Since some 
Other time law is governing its process, it would appear to 
Bevan interface controlled process. 

This, however, creates an anomaly in that mass transport 
was observed across this small strip of originally formed 
intermetallic which certainly had some sort of interface 
associated with it and yet the interfaces of the strips on 
each side of the region of mixing limit the intermetallic's 


growth to some small time power law. 
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A possible explanation for this behavior is that surface 
diffusion is the dominant mechanism. Instead of tellurium 
atoms moving across the interface, consider mass transport 
across the surface of this small strip of initially stoi- 
chiometric Bi2zTes;. This causes tellurium accumulation on 
the bismuth side and mass depletion on the tellurium rich 
side of the strip. Followed by diffusion into the film 
perpendicular to the substrate, the stoichiometric strip 
Will grow laterally. However, as the strip grows in lateral 
extent, the distance of aeitweten surface transport would 
increase. If the rate of transport across the surface 
decreases as the extent of the strip grows, as is experi- 
mentally observed when the extent of tellurium transport is 
limited to the region of mixing, then the rate of formation 
of the intermetallic is fast initially falling at an unknown / 
decreasing rate. 

This strip as it grows would also provide a barrier to 
the tellurium moving away from the tellurium rich side of 
the region mixing, which at the same time, through telluriun 
depletion, is creating strips of intermetallic, even though 
very thin in thickness normal to the substrate. 

With the intermetallic extending throughout the region 
of mixing, and now due to its extent in lateral dimension, 
and its respective interfaces, we would not expect and do 


not observe any growth with respect to time. 
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In summary, we can conclude that thin film couples of 
bismuth and tellurium form the intermetallic compound in the 
region of mixing and further growth of the intermetallic 
compound was not observed. Therefore, the intermetallic 
compound acts as a diffusion barrier and transport of bismuth 
mee tellurium through the intermetallic is interface con- 
trolled following some small time power law which was not 
decermined in this work. Due to the vapor pressure of 
tellurium, it is recommended that any further studies in 
the bismuth-tellurium system be done in a high vacuum condi- 


tion and/or temperatures below 200°C. 
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APPENDIX B 
EXPERIMENTAL RELATIONSHIPS USED FOR DEPOSITION OF 


ALLOY FILMS AND THIN FILM COUPLES 


Figure 38 shows the orientation relationship between 
the source, substrate and crystal monitor. The equations 
for evaporation from a small area (source) onto a parallel 
plane receiver had been previously reported [9]. The change 
in frequency (frequency shift) of the 5 MC crystal is related 
to the mass deposited on the crystal by equation (25) of the 
present work. The 8 variable represents the correction paranm- 
eter that adjusts for the orientation relationships of 
Figure 38. Therefore, the t in equation (25) represents the 


thickness of the film. 8 consists of the following variables 
Z 


Bees 1 
P= 7% Fe | cos20 oe) 


where Noma constant Of theverys tall 


R, = the substrate to source distance 
Ro = the crystal to source distance 
© = angle between the source normal and the crystal 


monitor. 


The thickness variation across the film is 


Coe i 
a Sa (31) 





he of 5 MC 
: ! CRYSTAL 
a See 
SUBSTRATE 
RAZOR BLADE R 
| R 
2 
STAINLESS @ 
STEEL SHEET 
SOURCE 
Rabie oes! Bye: 


Orientation Relationships between the Source, 


Substrate and Crystal Monitor 





Jest 


mere t = thickness of film directly above source 
Ce ecieChiessOt MimMieatnaistance GF from point on 
substrate directly above source. 

Using equation (31) the maximum thickness variation for a 
1,500 A film USdmne sthe values of r - 2 cm, R = 18 em, used 
meche present study, is 38 A, WilChabow Wied tere CLhrOv ron 
thickness measurement of the Dektac. 

Alloy deposition was carried out by 1) assuming that 
the film density was that of the reported bulk value and 
2) that the total thickness of the alloy film, t, was equal 
to the sum of the thickness of the two pure films. There- 


Here, let 


eo aes (32) 
Let desired composition of film be X weight% Ag, (1-X) 
weight% Au. 
X wel ht% Ag = _mass Ago x 100 
ee mass Ag + mass Au 
Pee Nout 
Pag tag * * Pau tau 4 Cs 


where A = substrate area on which gold and silver is 


deposited. 


Substituting equation (32) into equation (33) and solving 


or ty 
: tp, ("7% Ag) 
tag ~ Da (™® Ag) * 0,78 Aa) a 
g Py, (WS Ag Ag 
ancl fe =t-t (35) 


Au — Ag 





APPENDIX C 
COMPUTER PROGRAM USED IN ANALYZING DATA 


FOR THIN FILM STANDARDS 
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APPENDIX D 
COMPUTER PLOTTING PROGRAM WHICH PLOTS CONCENTRATION 
VERSUS RELATIVE DISTANCE FROM 


MICROPROBE INTENSITY DATA OF THIN FILM COUPLES 
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APPENDIX E 
COMPUTER PROGRAM WHICH NUMERICALLY SOLVES EXPERIMENTAL DATA 
OF THE INITIAL CONCENTRATION PROFILE FOR GIVEN VALUE 
OF D AND TIME 
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OTATISTICAL RELATIONSHIPS USED IN THE PRESENT WORK 
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APPENDIX F 


STATISTICAL RELATIONSHIPS USED IN THE PRESENT WORK 


Consider curve of the type Y = a,;X. Assume that for 
the data point (Yi, Xi) each Yi is an observation from a 
normal population with a mean a, Xi and constant variance 
o*. Assuming that measurement of Y locates it within a 
small interval of extent AY, the probability of observing 


the value Yi is from the normal distribution: 
N 
Sua WP (Yi Saya) - (36) 
i=l 


For S to be a minimun dome 
da, 





95 _ = - 25YX + 2a,EX2 = 0 
da, 
Therefore, 
eee PO! 
ai. yX2 573) 


The variance about the regression line is estimated by the 


square of the standard error of the estimate 


zs 4 2. 2 
SZ, ae OCC ae 3 oe (38) 


where N = the number of data points. 
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Lines drawn parallel to the regression line at vertical 
distance (Y) of + Syx> + 2Syy and + 3Syy will enclose 68%, 
95% and 99% of the data points. The value of Syy gives a 


measure of the dispersion about the mean. 


The coefficient of correlation is defined as 


Seale 
r = F a (39) 





where 


iMecwecOcthictenteot Correlation iS a confidence test on 
Biewlinearity of the data. For example, if the coefficient 
of correlation equals 0.658 and 0.708 we can say with 98 and 
Soo certainty that there is linear correlation of the data. 

If we consider a curve of the type Y = apg + a,X, then 


following the procedures above we can derive the following 


relationships 
L TOReR ORE (ORG OR 88 
* NEX2 - (EX)2 va 
Bm ND Xoo ane) 
DS Ser omGRC a a) 
Ss the square of the standard error of the estimate equals 


YX? 


ee INS SEO Siok 
me Gee (42) 
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iitemvariance of slope, V = Gh Ors (43) 
is 


The confidence band of the slope, say equal to b, is a fan 
Shaped area converging on the mean value of X and Y with 

slopes b + Vt, where t is the probability level determined 
from an F test on the coefficient of correlation (equation 


ao). 
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